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Q ■ Abstract 

(N : 

A simple and algorithmic description of matrix shape invariant potentials is pre- 
sented. The complete lists of generic matrix superpotentials of dimension 2x2 and 
of special superpotentials of dimension 3x3 are given explicitly. In addition, a con- 
structive description of superpotentials realized by matrices of arbitrary dimension is 
presented. In this way an extended class of integrable systems of coupled Schrodinger 
fH ■ equation is classified. Examples of such systems are considered in detail. New in- 

Qh' tegrable multidimensional models which are reduced to shape invariant systems via 

rl ! separation of variables are presented also. 
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1 Introduction 

Supersymmetry presents effective and elegant tools to solve quantum mechanical problems 
described by integrable Schrodinger equations. Unfortunately the class of known problems 
which can be solved using their supersymmetry is rather restricted since they should have 
the additional quality called shape invariance [1], and this feature appears to be rather 
rare. The classification of shape invariant (scalar) potentials is believed to be completed 
at least in the case when they include an additive variable parameter j2]. 

However there exist an important class of shape invariant potentials which is not classi- 
fied yet, and they are matrix valued potentials. Such potentials appear naturally in models 
using systems of Schrodinger-Pauli equations. A famous example of such model was pro- 
posed by Pron'ko and Stroganov (PS) [3J, its supersymmetric aspects were discovered in 
papers |1] and [5] . We note that there exist a relativistic version of the PS problem which 
is shape invariant too [B]. 

Examples of matrix superpotentials including shape invariant ones were discussed in 
[Z], |8], [9] [lO], [11]. A rather general approach to matrix superpotentials was proposed 
in paper [12], which, however, was restricted to their linear dependence on the variable 
parameter. 

A systematic study of matrix superpotentials was started in recent paper |T3] where we 
presented the complete description of a special class of irreducible matrix superpotentials. 
These superpotentials include terms linear and inverse in variable parameter, moreover, 
the linear terms where supposed to be proportional to the unit matrix. In this way we 
formulated five problems for systems of Schrodinger equations which are exactly solvable 
thanks to their shape invariance. Three of these problems are shape invariant with respect 
of shifts of two parameters, i.e., posses the dual shape invariance [15] . 

The present paper is a continuation and in some sense the completion of the previous 
one. We classify all irreducible matrix superpotentials realized by matrices of dimension 
2x2 with linear and inverse dependence on variable parameter. As a result we find 17 matrix 
potentials which are shape invariant and give rise to exactly solvable problems described by 
Schrodinger-Pauli equation. These potentials are defined up to sets of arbitrary parameters 
thus the number of non-equivalent integrable models which are presented here is rather 
large. They include as particular cases all superpotentials discussed in [5], [9], [12] and 
|13] . but also a number of new ones. Moreover, the list of found shape invariant potentials 
is complete, i.e., it includes all such potentials realized by 2 x 2 matrices. 

In addition, we present a constructive description of superpotentials realized by matrices 
of arbitrary dimension. The case of matrix superpotentials of dimension 3 x 3 is considered 
in more detail. A simple algorithm for construction of all non-equivalent 3x3 matrix 
superpotentials is presented. A certain subclass of such superpotentials is given explicitly. 

The found superpotentials give rise to one dimensional integrable models described by 
systems of coupled Schrodinger equations. However, some of these systems are nothing 
but reduced versions of multidimensional models, which appear as a result of the separa- 
tion of variables. Examples of such multidimensional systems are presented in section 8. 



These systems are integrable and most of them are new. In particular, we show that the 
superintegrable model for vector particles, proposed in |21], possesses supersymmetry with 
shape invariance, and so its solutions can be easily found using tools of SUSY quantum 
mechanics. The same is true for the arbitrary spin models considered in paper [21] , but we 
do not discuss them here. 

We also analyze five of found integrable systems in detail and calculate their spectrum 
and the related eigenfunctions. In particular we give new examples of matrix oscillator 
models. 

The paper is organized as follows. In section [2] we discuss restrictions imposed on 
superpotentials by the shape invariance condition and present the determining equations 
which should be solved to classify these superpotentials. In sections [3] and S] the 2x2 
matrix superpotentials are described and the complete list of them is presented. The 
case of arbitrary dimensional matrices is considered in section |5l the list of 3 x 3 matrix 
superpotentials can be found in section [61 

Sections [7| and [8] are devoted to discussion of new integrable systems of Schrodinger 
equations which has been effectively classified in the previous sections. In addition, in 
section [H] we discuss SUSY aspects of superintegrable models for arbitrary spin s proposed 
in ED. 



2 Shape invariance condition 

Let us consider a Schrodinger-Pauli type equation 

H^ij = E^iIj (2.1) 

where 

Hk = -^,+V,{x\ (2.2) 

and Vk{x) is a matrix- valued potential depending on variable x and parameter k. 

We suppose that Vk{x) is an n x n dimensional hermitian matrix, and that Hamiltonian 
Hk admits a factorization 

H, = ata- + c, (2.3) 

where 

«; = !: + w., at = -A + w., 

Ck is a constant and Wk{x) is a superpotential. 

Let us search for superpotentials which generate shape invariant potentials Vk{x). It 
means that Wk should satisfy the following condition 

w', + w;, = wl^^ - wUa + c, (2.4) 



were Ck and a are constants. 

In the following sections we classify shape invariant superpotentials, i.e., find matrices 
Wk depending on of x,k and satisfying conditions (12. 4p . More exactly, we find indecom- 
posable hermitian matrices whose dependence on k is is specified by terms proportional to 
k and p 

Let us consider superpotentials of the following generic form 

Wk = kQ + ^R + P (2.5) 

where P, R and Q are n x n Hermitian matrices depending on x. 

Superpotentials of generic form (12. 5p were discussed in paper [13] were we considered 
the case of arbitrary dimension matrices but restrict ourselves to the case when Q = Q{x) 
is proportional to the unit matrix. Rather surprisingly this supposition enables to make 
a completed classification of superpotentials (12. 5p . All such (irreducible) superpotentials 
include known scalar potentials listed in pj and five 2x2 matrix superpotentials found in 



To complete the classification presented in [13j let us consider generic superpotentials 
(12. 5p with arbitrary hermitian matrices Q, P and R. We suppose Wk be irreducible, i.e., 
matrices R, P and Q cannot be simultaneously transformed to a block diagonal form. 

Substituting (12. 5p into (12.41) . multiplying the obtained expression by k'^{k + a)^ and 
equating the multipliers for same powers of k we obtain the following determining equations: 

Q' = a{Q^ + uI), (2.6) 

P'--{Q,P} + >^I = 0, (2.7) 

{R,P} + XI = (2.8) 

R^ = u^I (2.9) 

where Q = ^Q, Q' = -^, {Q, P} = QP + PQ is an anticommutator of matrices Q and 
P, I is the unit matrix and x. A, u are constants. 

Equations (I2.6p - (l2.9p have been deduced in [13] where the anticommutator {Q, P} was 
reduced to doubled product of Q with P since Q was considered to be proportional to the 
unit matrix and so be commuting with P. 

The system (I2.6p - (l2.9p for generic matrices Q, P and R is much more complicated then 
in the case of diagonal Q. However it is possible to find its exact solutions for matrices of 
arbitrary dimension. 

3 Determining equations for 2x2 matrix superpoten- 
tials 

At the first step we restrict ourselves to the complete description of superpotentials (12. 5p 
which are matrices of dimension 2x2. In this case it is convenient to represent Q as a 



linear combination of Pauli matrices 

Q = goo-Q + qiCTi + g20"2 + gsCTs (3.1) 

where ctq = / is the unit matrix, 

ai=(j J), cr,= (j -^y a3=(j _° J . (3.2) 

Let is show that up to the unitary transformation reahzed by a constant matrix the 
matrix Q can be transformed to a diagonal form. Substituting (13.1 p into f l2.6p and equating 
coefficients for the linearly independent Pauli matrices aa, a = 1, 2, 3 we obtain the following 
system: 

q'^ = 2aqoqa, a = 1,2, 3, (3.3) 

It follows from (13. 3p that 

Qa = CaF{x), F{x) = cxp ( 2a qodx j (3.4) 

where c^ are integration constants. Since all g^ are expressed via the same functions of x 
multiplied by constants, we can transform Q to the diagonal form: 

rt _ / *?+ 



Q ^ UQU^ = (^ ^^ ^_ j (3.5) 

where g± = go ± cF{x), c = A/cf+cf+c|, and U is the constant matrix: 
c + cs - i(j2Ci + iaiC2 



U 



2^C{C + C3) 



In accordance with f l3.5p equation fl2.6p is reduced to the decoupled system of Riccati 
equations for q±: 

q'^ = a{ql + z/) (3.6) 

which is easily integrable. The corresponding matrices P can be found from equation (12. 7p : 

with p± being solutions of the following equation 

p'± = ap±q± + X, (3.8) 



and 

p = fxexp(-a l{q+ + q-)dx] (3.9) 

where /i is an integration constant and the asterisk denotes the complex conjugation. More- 
over, up to unitary transformations reahzed by matrices commuting with Q (13. 5p the con- 
stant /i can chosen be real and so we can restrict ourselves to p satisfying p* = p. 

Consider the remaining equations (12. 8 p and (12. 9p . In accordance with (12. 9 p R should 
be a constant matrix whose eigenvalues are ±u. Thus it can be represented as R = 
riO"i + r2cr2 + r^a-^ where ri, r2 and r^ are constants satisfying r\ + r^ + r\ = u"^, or 
alternatively, R = ±ul. Let w 7^ then, in order equation (12. 8p to be satisfied we have 
to exclude the second possibility and to set ri = h = p± = 0. As a result we obtain the 
general solution of the determining equations (I2.6l) - (l2.9p with w 7^ in the following form: 

P = aip, R = r^as + r2cr2, Q = q+a+ + q-a_ (3.10) 

where cr-t = (1 ± o"3)/2, q± are solutions of Riccati equation (13. 6p . p is the function defined 
by (13.91) and r^ are constants satisfying r^ + r^ = uj^ . 

If w = then conditions (12. 8p and (12.91) became trivial. The corresponding matrices Q 
and P are given by equations (13. 5p and (13.71) . 

4 Complete list of 2 x 2 matrix superpotentials 

Let us write the found matrix superpotentials explicitly and find the corresponding shape 
invariant potentials. 

All nonequivalent solutions of equations (13.61) are enumerated in the following formulae: 

g. = 0, z/ = 0, (4.1) 

(la = , 1^ = 0, 

A A^ 

g^ = -tan(A2; + c<^), z/ = — > 0, (4.2) 

a a^ 

A A^ 

g<j = tanh(Ax + Co-), v = < 0, 

a a^ 

A A^ 
qa = coth(Aa; + c^), v = < 0, 

A A^ ^'-'^ 

qa = , I' = ^ < 

where a = ±,c„ = ±c and c is an integration constant. 

Going over solutions (I4.ip - (l4.3p corresponding to the same values of parameter u it is 
not difficult to find the related entries of matrix P (13. 7p defined by equations (13. 8 p and 



( I3.9p . As a result we obtain the following list of superpotentials: 

WJ^^^ = A (k (cr+ tan(Ax + c) + (T_ tan(Aa; — c)) 



(4.4) 



+/icriA/sec(Ax — c) sec(Ax + c) -\ — R 
iy^2) ^ _)^ {-K{a+ coth(Aa; + c) + a. coth(Ax - c)) 

+/icri v^csch(Ax — c) csch(Ax + c) -\ — R J , 
W^^'^ = A (-/€(o"+ tanh(Ax + c) + a^ tanh(Ax - c)) 

+/icri ^sech(Aa; — c) secli(Ax + c) H — i? ) , 

W^i^^ = A (-fi:(o-+ tanh(Ax + c) + (7+ coth(Ax - c)) 

+/icri A/sech(Aa; + c) csch(Ax — c) H — R ] , 

K J 

W^^^ = A ( — fi;(cr+ tanh(Ax) + cr_) + /i(TiA/sech(Ax) exp(— Aa;) H — i? j 

Vl^^6) ^ _)^ ( -fi;(cr+ coth(Ax) + cr_) + /iCTi A/cscli(Aa;) exp(-Ax) + -i?J , (4.9) 

WP = -. (^^ + -^ + -^^L= + U (4.10) 

\x + c x-cj ^x^ - c2 K ^ ^ 

Wf^ = -K— + /iai^ + -i? (4.11) 

X \/X K 



(4.5) 



(4.6) 



(4.7) 



(4.^ 



PFf ) = A (-kI + /iexp(-Ax)ai - -ag) . (4.12) 

Here a± = |(o"o ± a^), R is the numeric matrix given by equation (I3.10p . k, /i and A are 
arbitrary parameters. 

Formulae (I4.4I) - (I4.11I) give the complete list of superpotentials corresponding to non- 
trivial matrices R. In particular this list includes superpotentials with Q being proportional 
to the unit matrix which has been discussed in paper [13]. These cases are specified by 
equation ( KT2\i and equations dOD, dSD, O, flCTI]) with c = and R = uai. 

Finally, let us add the list (I4.4p - (l4.12p by superpotentials with R = 0. Using equations 
(13. 5p . (I4.ip -( l4l3|) and (13. 7p . (13. 9 p we obtain the following expressions for operators (12. 5p : 

W^^^'^ = A (a-+(Ktan(Ax + c) + z/sec(Aa; + c)) 

\ (4.13) 

+cr_(fi;tan(Ax — c) + r sec(Ax — c)) + ficri^/sec{Xx — c) sec(Aa; + c) 1 , 

Wj:^^^ = -A ((T+(/tcoth(Ax + c) + z/csch(Ax + c)) 

/ r\ (4-14) 

+0"- (k coth(Ax — c) + r csch(Ax — c)) + /icri y csch(Ax — c) csch(Ax + c] 



iy^i2) ^ _^ {a+{Ktanh{Xx + c) + iysech{Xx + c)) 

+a_ (k coth(Aa; — c) + t csch(Aa; — c)) + /iO"! A/sech(Aa; — c) csch(Aa; + c) j , 
W^^'-^^ = -A (a-+(/ttanh(Ax + c) + z/sech(Ax + c)) 



(4.15) 



/ \ (4.16) 

+o"_(Ktanh(Ax — c) + r sech(Ax — c)) + /i(7iA/sech(Ax — c) sech(Ax + c) j , 

VF^^"^^ = —A fcr+(/ttanh Ax + z/sech Ax) + cr_K) + yUcriA/sech Axexp(— Ax) j , (4.17) 

WJ.^^' = —A (/t((T+(fi;coth Ax + z/csch Ax + (T_/t) + fiai A/csch Ax exp(— Ax) j , (4.18) 

H'™ = --. (^ + ^(x + c)) -a(^-^^ ^(x - c)) + ^42=. (4.19) 
Vx + c2 y Va;-c2 / y/x^ - c^ 

...(•17-1 /2fi;+l wxA /OJX \ 1 , , 

'*'' ' = -"- [-2r - TJ + ^- (t + ') - ''^■^- '^-2°) 

Formulae (I4.4p - O4.20p give the complete description of matrix superpotentials realized 
by matrices of dimension 2x2. These superpotentials are defined up to translations 
x—J-x + c, «;—)■«; + 7, and up to equivalence transformations realized by unitary matrices. 
In (I4.4p - (l4.20p we introduce the rescalled parameter k = - such that the transformations 
fc — !■ A;' = A; + a is reduced to: 

K ^ k' = K + 1. (4.21) 

The list (I4.4p - (l4.20p includes all superpotentials obtained earlier in [12j and [13], but 
also a number of new matrix superpotentials. The corresponding shape invariant potentials 
are easily calculated starting with superpotentials (14.40 - 04.200 and using the following 
definition: 

^» = iy«2_^«'^ , = 1,2,. ..,17. (4.22) 

To save a room we will not present all potentials (14.221) explicitly but restrict ourselves 
to discussions of particular examples of them, see section [71 

5 Matrix superpotentials of arbitrary dimension 

Let us consider generic superpotentials (12.51) with arbitrary dimension hermitian matrices 
Q, P and R. In this case we again come to the determining equations (I2.6p - (I2.9I) where 
Q, P and R are now hermitian matrices of dimension K x K with arbitrary integer K. 

In accordance with (12. 9p i? is a constant matrix whose eigenvalues are ±u). Thus up to 
unitary equivalence it can be chosen in the form 

R = uj(^^ ^ ), n + m = K (5.1) 



where /„ and Im are the unit matrices of dimension n x n and m x m correspondingly, n 
and m are the numbers of positive and negative eigenvalues of R. 
It is convenient to represent matrix Q in a block form: 

« = ( B. C ) (=-2) 

where A, B and C are matrices of dimension nxn, nxm and mxm. Using the analogous 
representation for P and taking into account relations (I2.8p we right it as 

where r = — tt"- 

Substituting fl5.1l) - fl5.3p into (12.61) and (12.71) we obtain the following equations for the 
block matrices: 

A' = a{A^ + BB^ + uQ, (5.4) 

C' = a{C^ + B^B + uIm), (5.5) 

B' = a{AB + BC), (5.6) 

P' = ^{AP + PC), (5.7) 

2tA + BP^ + PB^ = 2714, (5.8) 

~2tC + B^P + P^B = 2-JiIm (5.9) 

where 11= -. 

Thus the problem of description of matrix valued superpotentials which generate shape 
invariant potentials is reduced to finding the general solution of equations (I5.4p - (I5.9I) for 
irreducible sets of square matrices A, C and rectangular matrices B and P. Moreover, A 
and C are hermitian matrices whose dimension is n x n and mxm respectively while 
dimension of B and P is {nxm). Without loss of generality we suppose that n < m. 

The system (I5.4p - (l5.9p is rather complicated, nevertheless it can be solved explicitly. To 
save a room we shell not present its cumbersome general solution but restrict ourselves to 
the special subclass of solutions with trivial matrices B. In this case Jl = t = (otherwise 
the corresponding superpotentials are reduced to a direct sum of 2 x 2 matrices considered 
in the above), and the system is reduced to the following equations: 

A' = a{A^ + uQ, C = a{C^ + uIm), (5.10) 

P' = ^{AP + PC). (5.11) 

Without loss of generality the hermitian matrices A and C which solve equations (I5.10p 
can be chosen as diagonal ones, see Appendix. In other words their entries Aab and Cab 
can be represented as: 

Aab = SabQb, Cab = ^abQn+b (5.12) 



where 6ab is the Kronecker symbol and q„ [a = b or a = n + b) are solutions of the scalar 
Riccati equation 

g; = a{ql + u) (5.13) 

which is a direct consequence of fl5.10|) and fl5.12|) . Solutions of equations fl5.13|) are given 
by equations (|H])-(|43D. 

Thus matrices A, C and R are defined explicitly by relations fIS.lUI) . fl4.ip -( H^ and 
fIS.ip while matrices B are trivial. The remaining components of superpotentials f l2.5p are 
matrices P whose entries Pab are easy calculated integrating equations (15.71) : 

Pab = l^ab exp ( -a hqa + qn+b)dx j (5.14) 

where Hab are integration constants satisfying Hab = if^ba)*, and q^ with a = a,n + b are 
functions (l4.ip -f Hl3|l corresponding to the same value of parameter u/a. 

In analogous way we can describe a special subclass of matrix superpotentials (12. 5 p with 
trivial matrices R [14j. In this case it is convenient to start with diagonalization of matrix 
Q and write its entries as 

Qaa = SaaQa, a, a = 1,2, . . . , K (5.15) 

where q^ are functions satisfying equation (I5.13p . Then the corresponding entries of matrix 
P satisfying (15. 7p are defined as: 



Paa = f^acT exp ( -a [q^ + qa)dx j . (5.16) 

Functions g^ and q„ included into (I5.16P have to satisfy equation (I5.13P with the same value 
of parameter u. In addition, the matrix whose entries are integration constants fiaa should 
be hermitian. 



6 Superpotentials realized by 3x3 matrices 

Let us search for superpotentials (12.50 realized by matrices of dimension 3x3. We will 
restrict ourselves to the case when parameter u in the determining equation (12. 6p is equal 
to zero and find the complete list of the related superpotentials. Like (14.101) . (14.111) . (I4.19P 
and (I4.20p they are liner combinations of power functions of x + q with some constant q. 
There are three versions of the related matrices R whose general form is given in (15. ip : 

^-(t ;)' ^'-'^ 

R = h, (6.2) 

R = 03. (6.3) 



Let us start with the case presented in (16.11) . The corresponding matrices Q and P are 
given by formulae (15. 2p and (15. 3 p with 

where 01,02,03,0,61,62,^1 and p2 are unknown scalar functions. Moreover, 01,03 and c 
should be real, otherwise Q is not hermitian. Without loss of generality we suppose that 
Pi and P2 be imaginary, since applying a unitary transformation to Q, P and R these 
functions always can be reduced to a purely imaginary form. Moreover, the corresponding 
transformation matrix is diagonal. 

In the previous section we apriori restrict ourselves to trivial matrices B. Let us show 
that this restriction is not necessary at least for the considered case z/ = 0. 

First let us prove that system (I5.4p - (l5.9p is compatible iff 71 = r = 0, and it is true for 
all versions of matrix R enumerated in (I6.ip - (l6.3p . Calculating traces of matrices present 
in (15. 8 p and (15. 9p we obtain the following relation: 

r(Tr A + TrC) + 77(m - n) = (6.5) 

where ?ti = 2, 3, for versions (16.11) . (16.20 . (16.11) respectively, and n = 3 — m. 

Differentiating all terms in (16. 5p w.r.t. x and using equations (15.41) . (15. 5p we obtain: 

r(TrA^ + TrC^ + 2TrB^B + v{m + n)) = 0. (6.6) 

Three the first terms in brackets are positive defined and we supposed that u = 0. If 
r 7^ all terms in brackets should be zero. If the trace of the square of hermitian matrix 
is zero, this matrix is zero too, the same is true for matrix B^B. Thus for r 7^ matrix 
Q (15. 2p is trivial. To obtain non-trivial solutions we have to set r = 0, then from (16. 5p we 
obtain that /I = also. 

Substituting (16. 4p into (15.80 and (15.90 (remember that /I = r = 0) we obtain the 
following relations: 

Pibl-hp2 = 0, P2bl-b2Pi = 0, pi(6i-6*) = 0, p^{b2-b;) = 0. (6.7) 

In accordance with (16. 7p there are three qualitatively different possibilities: 

(o) : pi= P2 = 0, (b) : P162 = P261, bl = bi, b^ = 62 (6.8) 

and 

(c) : 61 = 62 = 0. (6.9) 

In the cases (16. 8p the corresponding superpotentials are reducible. Indeed, in the case 
(a) the only condition we need to satisfy is equation (12. 6p . But matrix Q can be diagonalized 
(see Appendix) and so the related superpotential can be reduced to the direct sum of three 
scalar potentials. 

10 



The only possibilities to realize the case (b), which differs from cases (a) and (c) is to 
suppose that 

Pi = ap2 and 61 = ab2 (6.10) 

or pi = ahi and p2 = «&2 where a is a constant parameter. The second possibility is 
excluded since Pa can be proportional to ba only in the case when these functions are 
reduced to constants (otherwise Pa and b^ are linearly independent, compare equations 
(15. 6p and (15. 7p ). and so this possibility is reduced to (I6.10p also. 

But if conditions (I6.10p are realized the corresponding superpotential is reducible too 
since the transformation W — )■ UWW with 



U 



makes it block diagonal. 

Thus to obtain an irreducible superpotential (12.51) we should impose the condition (16. 9p . 
and our problem is reduced to solving the system (I5.10p . (15.111) where u = 0. Like in section 
[3] the 2x2 matrix A can be chosen diagonal, i.e., we can set 02 = in (16. 4p while the 
remaining (diagonal) entries of matrix Q can be denoted as ai = qi,a2 = q2,C = gs, 
compare with (I5.12p . In accordance with (I5.13P with z/ = 0, functions qi can independently 
take the following values: 






1 


or 
or 
or 


Qi = 
q2 = 

q^ = 


= 
= 
= 


gi — - 


X + Ci 

1 


52 — • 
Q3 = - 


X + C2 
1 



X + C3 



(6.11) 



where Ci,C2 and C3 are integration constants. The corresponding values of pi and p2 are 
easy calculated using equation (15.141) . As a result we obtain the following irreducible 
superpotentials: 



(6.12) 



W = {Si - 1)^^ + (Si - 1)^^ + (Si - 1)- 

^ ^ X + Ci ^ ^ X + C2 X 

+ S, , ^' + S2 , ^' + -{2Sl - 1), 

^Jx{x + Ci) a/x(x + C2) 1^ 

W = iS', - 1)^ + {Si - 1)^ + S,^ + ^2^% + -(2S| - 1), (6.13) 

X X + C y/x V X + C K ^ ^ 

W = (5? - 1)^ + {Si - 1)^ + S^J^ + 53^#= + -(2^1 - 1), (6.14) 

X + C X y/X ^x{x + c) /« ^ ' 

W = {Si -1)^ + S,c + S2-^ + -{2Sl - 1) (6.15) 

X y/X K ^ ' 

11 



where c, Ci, C2,/i,/ii and fj,2 are integration constants, and 

/O 0\ /O i\ /O 

5i = -i , ^2=0 0, 53 = i 01 (6.16) 






-i 


i 












are matrices of spin s = 1. 

Formulae fl6.12l) - fl6.15l) give the complete list of 3 x 3 matrix superpotentials including 
matrix R in form (16. ip . If this matrix is proportional to the unit one (i.e., if the version 
(16. 2 p is realized), the related matrix P should be trivial, see equation (12. 8 p with A = 
0. Diagonalizing the corresponding matrix Q we obtain the direct sum of three scalar 
potentials, i.e., the related superpotentials are reducible. 

In the case (16. 3 p we again can restrict ourselves to diagonal matrices Q whose entries 
are enumerated in (16. lip . The corresponding matrices P can be calculated using equation 
(I5.16p . As a result we obtain the following superpotentials: 



(6.17) 



W = (Si - 1)^^ + (Si - 1)^^ + (Si - 1)- 

^ ^ 'x + Ci ^ ^ X + C2 X 

, C ^1 , C ^2 , (^ /is 

^yx{x + cij ^yx[x + c2) \/{x + cij{x + C2j 

W = (Si - 1)- + {Si - 1)^^ + Si^L + 82^^^ + S-. ^' , (6 18) 

W = is', - 1)- + 5ic + ^3^^ + 32-^. (6.19) 

X \/ X \/ X 



Formulae (I6.12p - (l6.19p present the complete list of irreducible 3x3 matrix superpo- 
tentials corresponding to zero value of parameter z/ in (12. 6p . 

In full analogy with the above we can find superpotentials with u nonzero. In this 
case the list of solutions (16. lip is changed to solutions (I4.ip - (l4.3p with a = 1, 2, 3 and the 
same value of u for all values of a. The corresponding matrices P again is calculated using 
equation (I5.14p for non-trivial R and equation (12. 8 p for R trivial. 

7 Examples of integrable matrix models 

Thus we obtain a collection of integrable models with matrix potentials. The related 
superpotentials of dimension 2x2 are given by equations (I4.4p - (l4.20p . They are defined up 
to arbitrary constants c. A, /i, z/, . . . . In addition, in section [5] we present the infinite number 
of superpotentials realized by matrices of arbitrary dimension. So we have a rather large 
database of shape invariant models, whose potentials have the form indicated in (I4.22p . 

Of course it is impossible to present a consistent analysis of all found models in one 
paper. But we can discuss at least some of them. In this and the following sections we 
consider some particular examples of found models. 
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7.1 Matrix Hamiltonians with Hydrogen atom spectra 

Let us start with the superpotential given by equation (14.101) . In addition to variable 
parameter k it includes an arbitrary parameter /i, and two additional parameters, r2 and 
r^, which define matrix R f l3.10p . Moreover, /i^ H-rg 7^ 0, otherwise operator f l4.10p reduces 
to a direct sum of two scalar superpotentials. 

The simplest version of the considered superpotential corresponds to the case /i = 
and r^ = 0,r2 = to. Then with the unitary transformation Wk — > W^ = UWk U\ U = 
(1 + %az)j^2 we transform Wn to the following (real) form: 

^^._^_^. (7,1) 

X K 

The corresponding potential (I4.22p looks as: 

K = ^ — T^^ + ^ + ^- C^-2 

Ju r\j 

Shape invariance of potential (17. 2p is almost evident. Calculating its superpartner 

V+ = Wl + W'^ we easily find that 

V^ = K+i + C, (7.3) 

where 

Using this property (which makes our model be in some aspects similar to the non- 
relativistic Hydrogen atom) we immediately find spectrum of Hamiltonian (12. 2p with po- 
tential (17:^ : 

En = -^, (7.4) 

where N = k, + n, n = 0,1,2, .. . 

The ground state vector iPq{k,x) should solve the equation 

a-M^,x)=(— + wAMi^,x) = (7.5) 

where W,^ is 2 x 2 matrix superpotential (17.11) and iPq{k,x) is a two component function: 

^o(«:,x)=(^). (7.6) 

Substituting (17. ip and (17.61) into (17. 5p we obtain the following system: 

^' - -^ - -e = 0, (7.7) 

X K 
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<' -Lo^ = 0. (7.8) 

Substituting (p = -^' obtained from (17.81) into (17.71) we come to the second-order equa- 
tion for ^: 

r - -r - ^e = 0, 

whose normahzable solutions are: 



„^, /UX\ K,+ 1 ,„ , 

e = a;x^ir, (— j , ^^ = ^- (7.9) 

where Ki,{-) are modified Bessel functions. The first component of function (17. 6p is easy 
calculated using (17.81) : 



^ = KiK+ l)x^-^K, (^) - ujx-K,^^ (^) . (7.10) 

Solutions (17.61) . (17. 9p . (I7.10p are square integrable for any positive k and < x < oo. 
Solution ^„(/t;,x) corresponding to the rfi^ exited state can be obtained from the ground 
state vector using the following standard relation of SUSY quantum mechanics (see, e.g. 

my 

^„(fi;,x) = a+a++i---a++„_i^o (/« + «, a;). (7.11) 

It is not difficult to show that vectors (17. lip are square integrable too provided k is 
positive. 

In analogous manner it is possible to handle ground state wave functions corresponding 
to superpotential (14.101) with other values of parameters r2,r3 and /x. In general case (but 
for r2/i 7^ 0) these wave functions are expressed via products of exponentials exp (— ^), 
powers of x and Kummer functions Uau (^^)- We will not present the corresponding 
cumbersome formulae here. 

7.2 Matrix Hamiltonians with oscillator spectra 

Consider the next relatively simple model which corresponds to superpotential (I4.20p with 
c = 0. Denoting Wk = VF^ we obtain the following potential: 



Ak^ — 1 Lo'^x'^ A*^ / x^ 



K = W^ - < = a+ ——- + -— + ^ - (k + 1)- 



.2 „2 ,.X /o,.., .._ , .X .. (7-12) 



u) X fi uj\ (3ujx 2k, + 1\ fi 



^ ' 4 X 2 j\ A 2x J x/x 



It is easy to verify that the superpartner of V^, i.e., V^+ = W^ + W^ is equal to V^+i up 
to a constant term: 

V,^ = V^+i + LO. (7.13) 
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In other words, potential fl7.12p is shape invariant in accordance with the definition given 
in pp. Moreover, hke in the case of supersymmetric oscillator, the superpartner potential 
V^ differs from \4+i by the constant term cu which does not depend on variable parameter 
k. 

Using standard tools of SUSY quantum mechanics it is possible to find the spectrum of 
system (^^, ([22D with hamiltonian (TTT^ : 

En = nu, n = 0,1,2,... (7.14) 

which coincides with the spectrum of supersymmetric oscillator. 

The ground state vector is defined as a square integrable solution of equation fl7.5p . 
Substituting fl4.20p and f l7.6p into (17. 5p we obtain the following system: 

, /2k +1 UX\ 1^ ^ r. /^,-,^^ 



«'+f«-;7F--°- ("'') 



Changing in (I7:T5|1 . (JUS]) 

(^ = exp I — I (^, ^ = exp ( -^ 1 ^, (7.17) 



2 \ / 2 

UX \ _ ^ / UJX 



solving equation (17.160 for ^ and substituting the found expression into (I7.15P we obtain 
the second order equation for (p: 



_„ (UJX ^ n\ _, jji^ _ 



Its solutions are linear combinations of Heun biconfiuent functions: (p = Ci(pi + C2<^2 where 
(^1 = ifB(— a+, 0, a_, 6, ex), (p2 = x'^^^HB{ci+i^ia--,h^cx) (7-18) 

where 



a± = l±K, b= J. , c=- . (7.19) 

\/co 4 



Thus, in accordance with (I7.17p . (17.181) and (I7.15P we have two ground state solutions (17. 6p 
with 



ux^ 



(p = (pi = exp( — — ) HB{-a+, 0, a_, 6, ex). 



^ = ^i = exp( — —\ 1^ — /J^(-a+,0,a_,&,ex) (720) 



H — 7:— I Hsi—a+i 0, a_, 6, ex) 



2/i \ v^ 2 
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and 

ip = 1^2 = x''^^HB{a+, 0, a_, b, ex), 



I ^x'^\ I y2ujx ^^. . ^ , , 

t = to = exp Hr^[a<, U, a_, o, ex) 

^ 4 ; 1^ 4/x Bl +' ' ' ' ^ ^^21^ 

1 / 2fi; + 1 wxi , , 

H ;;— ^b{0'+, 0, a_, 0, ex) . 



2/i \ v^ 2 

Functions (17. 6p whose components are defined in (17.201) and (I7.2ip are square integrable 
for any real values of parameters k, yU and positive u. 

Notice that for integer n solutions (17.201) and (I7.2ip are linearly dependent. We will not 
present here the cumbersome expression of the second solution linearly independent with 
(I7.20p which can be easy found solving system (17.151) . (I7.16P for k integer. 

One more matrix superpotential generating the spectrum of supersymmetric oscillator 
is given by equation (14.191) . Setting for simplicity 5 = we obtain the corresponding 
Hamiltonian (I4.22|) in the following form: 

K = -ri^ + c ) + , , + {k + -]uj- aifix —■ -j + 



+ "("-'K"+(^ + "-(^ 

Like (I7.12p . potential (I7.22p satisfies the form-invariance condition written in the form 
(I7.13p . The spectrum of the corresponding Hamiltonian (12. 2p is given by equation (I7.14p . 
The ground state vectors, i.e., solutions of equation (17. 5p where W^ = Wk, is superpo- 
tential (14.191) with 6 = 0, are given by equation (17. 6p with components (f and ^ given 
below: 



if = ipi = exp (— 4a;(x + c)^) (c^ — x'^)'^{e — x)2Hc I a,b_, 6+, d,r; 



X + e 
2e 



e = ei = ^ exp {-Mx + ef) (c^ - x'r (^^c («, &- , h, d, r; ^^ (7.23) 



+ { K+ -j He ia,b-,b+,d,r;— — 

where Hc{. . . ) is the confluent Heun function, 

1 13 

a = -Auc^, b± = K±-, d = 2uc^, r = 26+c^w + 7:^."^ + ^ - /"^ (7.24) 

2 2 8 

There exist one more ground state vector for Hamiltonian (12. 2p with potential (I7.22p 
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whose components are 

1 / X ~\~ c 

ip2 = exp(-4u;(x + c)^) (c^ - x^)^{c- x)''Hc ia,-b^,b+,d,r;— — 

^2 = — exp (-Auj(x + cV) (c-xY (c(2cK-x)Hc I a,-b-,b+,d,r;- — -] (7.25) 
ex ^ ^ \ \ 2c J 

(X ~\~ c\ 
a,-b^,b+,d,r]— — j 



x^ — c^ 



where a,b±,d and r are parameters defined by equation fl7.24p . 

For u > and k > functions (17.231) and (I7.25P are square integrable on whole real 
axis. 

7.3 Potentials including hyperbolic functions 

An important model of ordinary (scalar) SUSY quantum mechanics is described by Schro- 
dinger equation with the hyperbolic Scarf potential 

K = -k{k - 1) sech2(x). (7.26) 

This model possesses a peculiar nature at integer values of the parameter k, namely, it is 
a refiectionless (non-periodic finite-gap) system which is isospectral with the free quantum 
mechanical particle. In addition, this model possesses a hidden (bosonized) nonlinear 
supersymmetry [13]. 

Let us consider shape invariant matrix potentials including (I7.26P as an entry. The 
corresponding superpotential can be chosen in the form (14. 8 p where /i = and C3 = 0, C2 = 
CO. Then with the unitary transformation W^ — ?■ VF^ = UW^ W^ , U = (l + ias)/V2 we 
transform it to the real form: 

W^ = X (-K{a+ tanh(Ax) + a_) + ^i-V (7.27) 

The corresponding potential looks as: 

K = AM -a+K{K - 1) sech^(x) - aia;(tanh(Aa;) + 1) + ^ + k^ J . (7.28) 

Potential (I7.28P satisfies the shape invariance condition (17. 3p with 

thus the discrete spectrum of the corresponding Hamiltonian (12. 2p is given by the following 
formula: 
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where 

n = 0,l,..., K + n<0, {K + n)'^>oo. (7.30) 

Conditions f l7.30p will be justified in what follows. 

To find the ground state vector we should solve equation (17.51) with W^ and iIjo{k,x) 
given by formulae (I7.27P and (17. 6p correspondingly. This equation is easy integrable and 
has the following normalizable solutions: 

cp = y -. (1 -y)^-l2Fi{a,b,c;y), 
Here 2-Fi(«, b, c; y) is the hypergeometric function, 

1 /— K LU 

c = 1 V K + uj , = c-\ 1 , a = b — K — 1, 

K 2 2k (^732) 

y = -(tanh Xx + 1). 

Wave functions for exited states can be found starting with (I7.32p and using equation 
(17. lip . In order to functions (I7.32p and the corresponding functions be square integrable, 
parameters k,u and n should satisfy condition (I7.30p . see discussion of normalizability of 
state vectors including the hypergeometric function in section 10 of paper [13] . 

Consider also superpotential (14.81) with a; = and fj, ^ 0: 

Wk = —A f K(cr+(tanh Ax + cr_) + ai/iA/sech Axexp(— Ax) j . (7.33) 

The corresponding potential (14.221) have the following form: 

Vk = A^ (— cr+/i;(fi; — 1) sech^ Ax + k^ 

2 1 ^ . ^ ^ , .. .^ Ax , 3 , \ (7.34) 



+/i sech Axexp(— Ax) + ain{2K — 1) exp — sechs Ax | . 

Solving equation (17. 5p with W^ and ipo{f^,x) given by formulae (I7.27P and (17. 6p we find 
components of the ground state vector: 



V = y ^^"''f^''\ ^^-a ^ = y%l-yr^2F,ia,b,c;y), 



-V/t2 + 2/i2, y = -(tanh Ax + 1) 



a = z/— — , b = a + K + -, c=l + 2u, v 



which are square integrable for k < 0. The discrete spectrum of Hamiltonian (12. 2p with 
potential (17. 2p is given by the following formula: E = — A^(k + n)"^ where n are natural 
numbers satisfying the condition k + n < 0. If this condition is violated, the related 
eigenvectors (17. lip are not normalizable. 
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8 Multidimensional integrable models 

The models considered in previous subsections are one dimensional in spatial variable. Of 
course it is more interesting to search for multidimensional (especially, three dimensional) 
models which can be reduced to integrable models by separation of variables. Famous ex- 
amples of such models are the (non-relativistic) Hydrogen atom and the Pron'ko-Stroganov 
system [3] which can be reduced to a scalar and matrix shape invariant systems correspond- 
ingly. A more "fresh" example is the reduction of the AdS/CFT holographythe model to 
the Poschl- Teller system proposed in [TB] . 

In this section we consider new examples of the three-dimensional Schrodinger-Pauli 
equations which can be reduced to a shape invariant form by separation of variables. 
Moreover, the related effective potentials in radial variable belong to the shape invariant 
potentials deduced above. 

8.1 Spinor models 

Consider shape invariant potential generated by the following superpotential: 

W = (/ias - J - 1)- + ^^^^^1. (8.1) 

X 2{j + 1) 

This operator belongs to the list of matrix superpotentials presented in section HI see 
equation (I4.10p . More exactly, to obtain (18. ip it is necessary to set c = 0, k = j + 1, r2 = 
and ri = — ^ in (I4.10p and (I3.10p . Then such specified superpotential Wk appears to be 
unitary equivalent to W (18. ip . namely, 

W = UWPU^ with U = ^(1 + ids). (8.2) 

V2 

Calculating potential (I4.22p corresponding to superpotential (18. ip with /i = | we obtain: 

V = W^-W'= (j(j + 1) + i - (^- + -3) ^ - ^a, (8.3) 

By construction, potential (18. 3p is shape invariant, thus the related eigenvalue problem 
92 



dx"^ 



+ V]'^ = EiP (8.4) 



can be solved exactly using standard tools of SUSY quantum mechanics. The corresponding 
ground state vector is a two-component function (17. 6p with 

^ = y^+lK,{y), ^ = y^+lKo{y), (8.5) 

where y = nrli) ■ Energy spectrum is given by equation (I7.4p with N = 2j + n + 1, n = 



0,1,... 



2{i+i) 
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The eigenvalue problem fl8.4p with potential f l8.3p includes the only independent variable 
X. However, it can be treated as a radial equation corresponding to the following three 
dimensional Hamiltonian with the Pauli type potential: 

i7 = -A + wcr • B, B = ^. (8.6) 

x^ 

Here 5 is the Laplace operator, ^ cr is a spin vector whose components are Pauli matrices 
f l3.2p . and B is the coordinate three vector divided by x^ = x\ + X2 + x\. 

Of course, B has nothing to do with the magnetic field since V ■ B 7^ 0. However, it can 
represent another field, e.g., the axion one [17]. Existence of such solutions for equations 
of axion electrodynamics was indicated recently [T8] . 

Expanding solutions of the eigenvalue problem for Hamiltonian f l8.6p via spherical 
spinors we obtain exactly equation (18. 5p for radial functions. We will not present the 
corresponding calculation here which can be done using the standard representations for 
the Laplace operator and matrix <t ■ x in the spherical spinor basis, which can be found, 
e.g., in [E]. 

Let us remind that the Pron'ko-Stroganov model is based on the following (rescalled) 
Hamiltonian 

if = -A + ^i^^-^, r' = xl + xl (8.7) 

which is reduced to the following form in cylindrical variables: 

92 1 1 

Hm = -^ + m{m-a:i)—^ai- (8.8) 

(we ignore derivatives w.r.t. X3). 

Hamiltonian (18. 8p can be expressed in the form (12. 3p . Moreover, the corresponding 
superpotential again can be obtained starting with superpotential (I4.1U|) by setting k = 
m+|,/i = |,c = r2 = 0,ri = | and making transformation (18. 2p . 

8.2 Vector models 

Hamiltonian (18. 7p corresponds to the physically realizable system, i.e., the neutral fermion 
moving in the field of straight line constant current. A natural desire to generalize this 
model for particles with spin higher then | appears to be hardly satisfied since if we simple 
change Pauli matrices in (18. 7p by matrices, say, of spin one, the resultant model will not 
be integrable [20] . 

In paper [2T] integrable generalizations of model Hamiltonian (18. 7p to the case of arbi- 
trary spin have been formulated. The price paid for this progress was the essential com- 
plication of the Pauli interaction term present in (18. 7p . However, there are rather strong 
physical arguments for such a complication [21], see also [22] for arguments obtained in 
frames of the relativistic approach. 
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In this section we present a new formulation of the spin-one Pron'ko model [2T] . Doing 
this we perform the following goals: to apply our abstract analysis of shape invariant matrix 
potentials to a physically relevant system and to show that this model is shape invariant 
and so can be easily solved using SUSY technique. 

Let us start with superpotential fl6.12p with ci = C2 = yU2 = 0, /ii = 1. Making the 
unitary (rotation) transformation W — )■ UWW with U = exp (15271/4) and changing the 
notations x ^ r, k, ^ m + ^ we reduce it to the following form: 



W = -S- 



1 ^ UJ /^ ^o ^\ 2m 



3 



- (25,^ - 1) - -^. (8.9) 



r 2m 

In addition, we transform the spin matrices to the Gelfand-Tsetlin form: 



^2 = -^ I 1 -1 I , ^3 = I I . (8.10) 







-1 





1 





-1 





1 






/I 











\0 


-1 


u' 





The corresponding shape-invariant potential looks as 

V^ = W^ + W'=[{m-S,r-\^y .{23^,-1)1 + j^£^^. (8.11) 

So far we simple represented one of the numerous supersymmetric toys classified in 
the above. Now we are ready to formulate a two dimensional model which generates the 
effective potential (18. lip . The corresponding Hamiltonian can be written as: 

d' d- 2(S-H)^-H^ 
^-'dxl-dxl \K\ ■ ^ ^ 

Here H is the two-dimensional vector of magnetic field generated by an infinite straight 
current; its components are Hi = q^ and H2 = ~q^- 

First we note that the last term in (I8.12p is a particular case of the interaction term 
found in [2l], see equations (15), (21), (29) here for s = 1 and /3i = /3o. However, we believe 
that our formulation of this term is more transparent. 

Secondly, introducing radial and angular variables such that xi = r cos 6, X2 = r sin 6, 
and expanding eigenfunctions of Hamiltonian (I8.12p via eigenfunctions ipm of the symmetry 

operator J3 = i [x2-£, ^i'£r ) "*" "^3 which can be written as: 

. /exp(i(m + 1)^)01 (r) 
V^^ = — exp(im^)02(r) | (8.13) 

\.exp(i(m — l)9)(f)^(r) 



,2 



we come to the following hamiltonian in radial variables: H = — ^ + Vm where Vm is the 
effective potential which coincides with (18. lip . Thus Hamiltonian (I8.12p is shape invariant 
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and its discrete spectrum and the corresponding eigenvectors are easily calculated using 
the standard tools of SUSY quantum mechanics. 

To end this section we present one more integrable model for vector boson. This model 
is three dimensional in spatial variables and is characterized by the following Hamiltonian 

H = -A + J^^-^l~^\ (8.14) 

|B| 

Here B is the three dimensional vector defined in (18.61) and S is the matrix vector whose 
components are given in equation (18.101) . 

Like (18. 6p . Hamiltonian (I8.14p corresponds to the shape invariant potential in radial 
variables, which looks as 

V = {jU + 1) + 5| - (2j + 1)^3) ^ + {2S', - 1)-. 

The corresponding superpotential can be obtained from (18. 9p changing m — )■ j ' + |. 

9 Discussion 

In spite of that (scalar) shape invariant potentials had been classified long time ago, there 
exist a great number of other such potentials which where not known till now, and they 
belong to the class of matrix potentials. The first attempt to classify these potentials which 
we made in recent paper [13] enabled to find five types of them which are defined up to 
arbitrary parameters. They give rise to new integrable systems of Schrodinger equations 
which can be easily solved within the standard technique of SUSY quantum mechanics. 

In the present paper we present an infinite number of such integrable systems. In 
particular we present the list of superpotentials realized by matrices of dimension 2x2, see 
equations (I4.4p - (l4.20p . The main value of the list is its completeness, i.e., it includes all 
superpotentials realized by 2 x 2 matrices which correspond to Schrodinger-Pauli systems 
(12. ip which are shape invariant w.r.t. shifts of variable parameters. 

In section[5]an extended class of arbitrary dimension matrix superpotentials is described. 
We do not present the proof that this class includes all irreducible matrix superpotentials. 
However such assumption looks rather plausible since by a consequent differentiation of 
equations (15.80 and (15. 9p with using conditions (I5.4p - (l5.7p we can obtain an infinite number 
of algebraic compatibility conditions for system (I5.4p - (l5.9p which are nontrivial but can 
be satisfied asking for B be equal to zero. An alternative solution of these compatibility 
conditions is P = 0, but it leads to reducible superpotentials. Computing experiments with 
system (I5.4p - (I5.9I) for the cases of n x ra matrix superpotential with ra < 5 also support the 
assumption 5 = 0, see section [6] where we did not make this assumption a priori but prove 
it. Notice that for n = 2 this condition is not necessarily satisfied, but it is seemed to be 
the only exceptional case. 

Nevertheless in section [5] we consider the condition i? = as an additional requirement, 
which enables to find superpotentials of arbitrary dimension in a straightforward way. 
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Thus we obtain an entire collection of integrable systems of Schrodinger equations. Some 
examples of these models are considered in section [7] where we we find their energy spectra 
and ground state solutions. Among them there are two oscillator-like matrix models whose 
spectra are linear in the main quantum number, see section 17.21 

One dimensional integrable models classified in the present paper are especially interest- 
ing in the cases when they can be used to construct solutions of two- and three-dimensional 
systems. A perfect example of such shape invariant system is the radial equation for the 
Hydrogen atom. Thus an important task is to search for multidimensional (in particu- 
lar, two- and three-dimensional) models which can be reduced to found shape invariant 
systems after separation of variables. Some results of our search can be found in section 
[S] were we present new integrable problems for two- and three-dimensional equations of 
Schrodinger-Pauli type. In particular we discuss SUSY aspects of the Pron'ko-Stroganov 
model generalized to the case of vector particles (such generalization was proposed in pa- 
per [21]). It happens that the spin-one model is shape invariant and so it can be easily 
integrated using tools of SUSY quantum mechanics. 

Except the case s = 1 we did not discuss superintegrable models proposed in [2T] for 
arbitrary spin s. However, it is possible to show that they are supersymmetric too. 

Let us note that the results presented in section [8] can be considered only as an adver- 
tisement, and we plane to present the detailed discussion of integrable multidimensional 
models in the following publications. 

A Solutions of matrix Riccati equations 

The corner stone of our classification of matrix superpotentials is the diagonalization of 
matrix Riccati equations (12.61) and (15.101) . For completeness, we present here a simple 
and constructive algorithm of such diagonalization for hermitian matrices of any dimension 
while the case of 2 x 2 matrices was already discussed in section [31 

Let us start with equation (12. 6p and make the following change of the dependent variable 
Q: 

Q = M + qI (Al) 

where I is the unit matrix and q = q^ is one of the particular solutions (14.21) of the 
scalar Riccati equation (I5.13P (the unnecessary subindex a can be omitted). Thanks to 
the randomness of integration constants in solutions (14.21) we always can suppose that the 
unknown matrix M be nondegenerated. 

Substituting (jXl]) into ([2SD and using the identity (M^^)' = -M-^M'M-^ we obtain 
the following /mear equation for M'^: 

(M~i)' = -a{I + 2qAr^). (A2) 

Its solutions have the following generic form: 

M-i = p{x)I + e{x)C (A3) 
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where p{x) and 6{x) are scalar real valued functions whose exact expressions can be easily 
found for any particular solution (14.21) . and C is an arbitrary constant matrix. 

In accordance with (]A3|) matrix M~^ is diagonalizable (remember that M~^ should be 
hermitian). The same is true for matrix M and so matrix Q (]A2|) is diagonalizable too. 

In complete analogy with the above one can justify the diagonalizability of matrices A 
and C satisfying equations (15.101) . 
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